Quadrupole Moments of Collective Structures up to Spin ~ 65H in 157 Er and 
158 Er: A Challenge for Understanding Triaxiality in Nuclei 
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Abstract 

The transition quadrupole moments, Q t , of four weakly populated collective bands up to spin ~ 65fi in 157 I58 Er have 
been measured to be ~11 eb demonstrating that these sequences are associated with large deformations. However, 
the data are inconsistent with calculated values from cranked Nilsson-Strutinsky calculations that predict the lowest 
energy triaxial shape to be associated with rotation about the short principal axis. The data appear to favor either 
a stable triaxial shape rotating about the intermediate axis or, alternatively, a triaxial shape with larger deformation 
rotating about the short axis. These new results challenge the present understanding of triaxiality in nuclei. 

Keywords: Ultrahigh-spin collective structures, transition quadrupole moments, cranked Nilsson-Strutinsky 

calculations, triaxial nuclear shape 
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The question of the occurrence of asymmetric shapes 
in nature is a general one touching on diverse sys- 
tems spanning many length scales, for example, large 
galaxies (diameter ~10 20 m), planets (~10 8 - 10 9 m), 
asteroids (~10 5 - 10 6 m), and small red blood cells 
(~10~ 5 m) or atomic clusters (~1(T 9 - 10~ 8 m). In nu- 
clear physics, the occurrence of stable asymmetric or 
triaxial shapes (~10~ 14 m) is a long-standing prediction 
of theory [1], and had been sought for experimentally 
for decades. However, it was not until 2001 that com- 
pelling evidence for nuclei with a robust triaxial shape 
was revealed in the rare-earth region (A~165) through 
the observation of rotational structures associated with 
the unique "wobbling" excitation mode [2]. A triaxial 
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nuclear shape, which has distinct short, intermediate, 
and long principal axes, is commonly described using 
the parameters (s2, y) of the Lund convention [ 3], where 
S2 and y represent the eccentricity from sphericity and 
triaxiality, respectively. At high spin, collective rotation 
about the short axis, corresponding to a positive y value 
(0° < y < 60°), usually has the lowest excitation energy 
based on moment of inertia considerations flU. Thus, 
this mode is expected to be favored over rotation about 
the intermediate axis (-60° < y < 0°). 

Recently, four rotational bands with high moments of 
inertia were observed at angular momenta up to 65ft in 
I57,i58g r Th ese sequences bypass the energetically 
favored "band-terminating" states [7, 8] and extend dis- 
crete y-ray spectroscopy to the so-called "ultrahigh-spin 
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regime" (/ = 50 - 70). Cranking calculations |0, [l^l 
had predicted collective configurations at large defor- 
mation to become competitive in energy for spins above 
5QH. Configurations with £2 ~ 0.34 and a positive 
value of y=20° - 25° were predicted to be low in en- 
ergy in Ref. [ 10] and re-confirmed in the calculations re- 
ported in Ref. Jg]. Thus, this shape minimum had been 
suggested to correspond to the ultrahigh-spin bands in 
I57,i58£ r This suggestion provided an opportunity to 
test the cranking model, a description widely employed 
in high-spin nuclear physics in the limit of extreme an- 
gular momentum. Hence, with the purpose of investi- 
gating the deformation associated with these structures 
at ultrahigh spin and of elucidating their character, a 
DSAM experiment was performed to extract transition 
quadrupole moments. 

In this letter, the transition quadrupole moments (Q t ) 
of the collective bands with high moments of inertia, 
bands 1 and 2 (adopting the labelling of Ref. in 
I57,i58j7 r ^ re p 0r i; ec i. In comparison with the calcula- 
tions mentioned above and new calculations using the 
cranked Nilsson-Strutinsky (CNS) model 111 ill , these re- 
sults appear to be inconsistent with the rotation of the 
theoretically favored triaxial shape. The implications of 
this rather surprising observation are discussed. 

The lifetime measurement was carried out using the 
Doppler Shift Attenuation Method (DSAM) 13. A 
215 MeV 48 Ca beam was delivered by the ATLAS fa- 
cility at Argonne National Laboratory and bombarded 
a 1 mg/cm 2 114 Cd target backed by a 13 mg/cm 2 197 Au 
layer. A 0.07 mg/cm 2 27 Al layer between Cd and Au 
was used to prevent the migration of the target material 
into the backing. The emitted y rays were detected by 
the Gammasphere spectrometer 11311 . which consisted 
of 101 Compton-suppressed HPGe detectors. A total 
of 9.9 x 10 9 events were accumulated, each one con- 
taining at least four coincident y rays. Due to the low 
intensity of the bands of interest, ~10~ 4 of the respec- 
tive channel intensity, several techniques 11141 Il5ll were 
utilized in the data analysis. An initial analysis using 
the RadWare software package II 1611 proved beneficial 
before sorting the data into a BLUE database II 1 711 to 
generate background-subtracted coincidence spectra at 
different detector angles lfl8ll . Examples of the resulting 
spectra at selected angles can be found in Fig.Q] 

The fractional Doppler shifts F(t) and the associ- 
ated errors for the inband transitions covering an esti- 
mated spin range of 30 - 6QH |6|] were subsequently ex- 
tracted from linear fits of the energy shifts as a function 
of detector angle 6. The resulting F(r) values for the 
four bands are presented in Fig. [2] The larger F(t) un- 
certainties of the two band 2 sequences relate to their 
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Figure 1 : Double-gated summed coincidence spectra from the present 
data for band 1 in l58 Er at two typical angles, 35° (top panel) and 145° 
(bottom) with respect to the beam direction. The triple-gated summed 
coincidence spectrum from the thin-target data |6] for the same band 
is shown as a reference in the middle panel. The energy values of 
the inband transitions (in keV) are given in the thin-target spectrum, 
while the peak positions of shifted y rays observed at the two angles 
are marked by lines in the other two spectra. 
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weaker intensities as compared with the bands 1 . The 
fit of the F(t) values was performed using the computer 
code MLIFETIME [19], which uses stopping powers 
provided by the SRIM 2003 code Q. It is worth men- 
tioning two of the assumptions made in the analysis: (1) 
all levels in a band have the same transition quadrupole 
moment, Q t ; and (2) the sidefeeding into each level in 
a band is modeled as a single cascade with a common, 
constant quadrupole moment, Q s f , and characterized by 
the same dynamic moment of inertia as the main band 
into which it feeds. Ax 2 minimization with the param- 
eters Q t and 2sf was then performed to the measured 
F(j) values for each of the four bands. This commonly 
adopted procedure has proved to be reliable for SD and 
TSD bands in this region, see Refs. 1121 1 12211 . for ex- 
ample. Details of the fitting process can be found in 
Refs. |l5, 23]. The resulting transition quadrupole mo- 
ments, Q t , and side-feeding quadrupole moments, Qsf, 
are summarized in the insets of Fig. [2] For compari- 
son, the ground-state band of 158 Er has a measured Q t 
of ~6 eb (for states up to spin 20ft) Il24ll and is asso- 
ciated with a normal prolate shape (S2 ~ 0.2, y=0°). 
Thus, these four bands with measured Q t ~l 1 eb are as- 
sociated with a large deformation. 

It is well known that an additional systematic error of 
~15% needs to be added for Q t measurements with the 
DSAM technique II25I1 . However, another measurement, 
focussing on I54 Er, was performed under very similar 
experimental conditions 12611. where the Q t of the yrast 
superdeformed band in 15 'Dy was extracted. The mea- 
sured value of 17+2 eb agrees with the reported one 
of 16.9*Q3 eb of Ref. 12111 and is reproduced well by 
the CNS calculations for a prolate superdeformed shape 
(see Fig. [_}. Hence, this comparison provides a con- 
sistency check for the experimental approach and the 
stopping powers used in the present data analysis. This 
151 Dy Q t value can also be seen in Fig. [3] where the data 
for the I58 Er bands are displayed and a comparison with 
CNS results is provided. 

In order to interpret these results, extensive CNS cal- 
culations without pairing [11] have been performed for 
I58 Er. The potential energy surfaces (PES) for states 
with every combination of parity and signature exhibit 
similar energy minima as a function of spin. PES plots 
for negative parity and a = 1 signature (odd spin) at 
two representative spin values of 49 and 69ft are given 
in Fig. Ufa) and SJb), respectively. In Fig. St a), there 
are four energy minima, one prolate, which is close to 
an axially symmetric shape with enhanced deformation 
at £2 ~ 0.28 (labelled as ED) when compared with the 
low-spin (7<30) yrast structure with E2 ~ 0.2 l24ll . a 
non-collective oblate minimum at £2 ~ 0.15 and y=60° 
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Figure 2: (Color online) Measured F(t) values as a function of the 
y-ray energy with best- fit curves for the four bands in 157 Er and 158 Er. 
The two horizontal dashed lines in each panel show the range of initial 
recoil velocities of the Er nuclei within the 114 Cd target layer. The 
insets summarize the measured quadrupole moments. The error bars 
are statistical only, i.e., they do not include the ~15% error associated 
with the systematic uncertainty in the stopping powers ['25] (see text). 
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Figure 3: (Color online) Measured transition quadrupole moments 
(Qt) of bands 1 and 2 in I58 Er and of the yrast SD band in 151 Dy (in- 
set), compared with the predicted Q t values associated with the min- 
ima of interest. The experimental results are plotted with diamonds 
for the present work and a square for Ref. |21]. The calculated CNS 
values in the spin range 30 - 60S are illustrated as horizontal shaded 
areas with the heights proportional to their ranges. See text for details. 
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Figure 4: (Color online) Potential energy surfaces in I58 Er calculated 
for negative parity and signature a = 1 at spins / = 49 (a) and 69 (b), 
respectively. The contour line separation is 0.25 MeV. The minima of 
interest are labelled (see text). 



(labelled as Yrast), and two triaxial ones with S2 ~ 0.34, 
y~20° (TSD1) and y~-20° (TSD2). In Fig. 12b), the 
higher-spin PES plot, another more deformed triaxial 
minimum at S2 ~ 0.43 and y~25° (TSD3) and a su- 
perdeformed one with £2 ~ 0.63 and y~0° (SD) emerge 
together with a higher deformation oblate yrast struc- 
ture. Of these minima, the non-collective oblate yrast 
minimum is associated with the terminating states Ollg], 
which are calculated lowest in energy in the spin range 
40 - 70/z (see Fig. 0. The positive-y triaxial minimum 
TSD1 is the lowest energy strongly deformed minimum 
throughout the spin range 30 - 657z. Therefore, it was 
the favored candidate for interpreting the ultrahigh-spin 
structures observed in Ref. 06J], prior to the present Q t 
measurement. 



The general formula of Q t for a triaxial nucleus [27] 
in terms of the density distribution along the three prin- 
cipal axes, x (rotation axis), y, and z (symmetry axis for 
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Figure 5: (Color online) Calculated excitation energies of the lowest 
states associated with the three triaxial-shape minima and of the yrast 
states in 158 Er plotted as a function of spin. The rotating liquid drop 
reference (Em) was calculated using the approach of Ref. [ll|]. 



y=0°), is 



Qt 



z 2 ) id, which is equivalent 



where Q 2 2(x) = iJUy 
to Eq. 15 in Ref. [28]. The same radius parameter, 
ro = 1.2 fm, has been employed as in previous stud- 
ies of axial SD nuclei, see Ref. I29I1 . for example. In the 
I58 Er case, the Q t values calculated with this formula 
are 6.8 - 7.9 eb for TSD1, 9.1 - 10.3 eb for TSD2, and 
8.6 - 9.9 eb for TSD3. These energy minima migrate 
gradually in the (s2, y) plane such that the associated Q t 
value for each minimum decreases slightly with increas- 
ing spin and forms a range, illustrated by a shaded area 
in Fig. [3] When calculating the Q t value for a triaxial 
shape, it has often been assumed that Q t depends on the 
y parameter as Q t = ie 2 o(z)|cos(y + 30°)/cos(30°) Q, 
where Q.2o(z) = (2z 2 - x 2 - y 2 ) yj] is the quadrupole 
moment of an axially deformed nucleus. The Q t val- 
ues for 158 Er calculated with this simple, but approxi- 
mate expression are smaller by ~8%, for the positive- 
y minima (TSD1 and TSD3), and larger by ~8%, for 
the negative-y minimum (TSD2), than those computed 
with the general formula above. However, these two 
formulae are equivalent for y=0°, and both give similar 
results for near-axially deformed minima, for example, 
g t ~6.7 - 7.9 eb for the ED case. 

The theoretical Q t values associated with all collec- 
tive minima, obtained with the general expression in the 
spin range 30 - 60H, are compared with the measured Q t 
values of the 158 Er bands 1 and 2 in Fig. [3] It is clear that 



the experimental Q t values, ~ 1 1 eb, are not consistent 
with either the calculated axially deformed minima, ED 
and SD, or the theoretically favored TSD1 minimum. 
The data are close to, but still larger, than the TSD2 
and TSD3 triaxial minima. Without a clear preference, 
TSD2 still has calculated Q t values slightly closer to the 
measured ones. However, TSD3 is predicted to be com- 
petitive with TSD2 in the 30 - 55ft spin range, and it is 
even lower than TSD2 for spins above 55ft (see Fig. O. 
Therefore, when comparisons of Q t are made with the 
CNS calculations, the ultrahigh-spin bands observed in 
158 Er appear most compatible with either the negative - 
y minimum TSD2 (s2 ~ 0.34) or the higher-deformed 
positive-y minimum TSD3 (S2 ~ 0.43). 

This is a surprising result, considering that the cal- 
culations consistently predict a 1 MeV gap between the 
TSD1 minimum and those associated with TSD2 and 
TSD3 in the 30 - 60ft spin range (Fig.|5]l- The observed 
bands appear to be the lowest in excitation energy since 
they are the only ones observed and the intensity pro- 
file indicates that they receive feeding from higher ly- 
ing states over the entire 40 - 60ft spin range. The in- 
formation on Q t values, predicted minima, and feeding 
patterns for the two 157 Er bands is very similar to that 
for 158 Er. Therefore, the same general conclusion can 
be drawn in that there is an inconsistency between the 
observations and the previously favored theoretical in- 
terpretation [6] for all the four bands observed in 157 Er 
and 158 Er. The lifetime data appear to rule out TSD1 
as the relevant minimum and, while the calculated Q t 
values for TSD2 and TSD3 are closer to the measured 
value for band 1, they are still outside the quoted errors. 

The present cranking calculations do not account for 
the Q t data satisfactorily, although it is noted that these 
calculations also carry some uncertainties. It would be 
interesting to see what other theoretical approaches pre- 
dict for the re-emergence of collectivity at these high 
spins in 157 Er and 158 Er. In the tilted-axis cranking 
(TAC) model lf3~lll . for example, it has been suggested 
that for a triaxial shape the rotational axis could tilt be- 
tween the intermediate and short principal axes. Thus, a 
rotational band may be associated with a mixing of the 
positive- and negative-y minima rf32ll and would have 
an average Q t value between the positive- and negative- 
y limits. Perhaps, the present observations should be 
linked with the differences of Q t values between the- 
ory and experiment for TSD bands and wobbling bands 
found in some Lu and Tm isotopes Ill5ll33ll34ll . While a 
satisfactory account can be given with CNS calculations 
for prolate superdeformed shapes ll29ll . the situation is 
different for bands which appear to be associated with 
triaxial strongly deformed shapes. On the other hand, 



the bands with large deformation in Hf have been de- 
scribed well i35ll in the CNS approach with a strongly 
deformed triaxial configuration. Obviously, a fully co- 
herent understanding is yet to emerge. 

In summary, DSAM measurements demonstrate that 
the ultrahigh-spin sequences in 157 158 Er are associated 
with strongly deformed shapes. The measured Q t val- 
ues appear to be more compatible, within the CNS the- 
oretical framework, with a negative-y triaxial deformed 
minimum (S2 ~ 0.34) or a positive-y minimum with 
larger deformation (£2~0.43) rather than with the en- 
ergetically favored positive-y triaxial shape (g2~0.34) 
predicted by theory. Discrepancies, however, still re- 
main between the experimental measurements and the 
theoretical predictions. This represents a challenge for 
the understanding of the triaxial degree of freedom in 
nuclei. 
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